Testing General Metric Theories of Gravity with Bursting Neutron Stars 
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I show that several observable properties of bursting neutron stars in metric theories of gravity 
can be calculated using only conservation laws, Killing symmetries, and the Einstein equivalence 
principle, without requiring the validity of the general relativistic field equations. I calculate, in 
particular, the gravitational redshift of a surface atomic line, the touchdown luminosity of a radius- 
expansion burst, which is believed to be equal to the Eddington critical luminosity, and the apparent 
surface area of a neutron star as measured during the cooling tails of bursts. I show that, for a general 
metric theory of gravity, the apparent surface area of a neutron star depends on the coordinate radius 
of the stellar surface and on its gravitational redshift in the exact same way as in general relativity. 
On the other hand, the Eddington critical luminosity depends also on an additional parameter that 
measures the degree to which the general relativistic field equations are satisfied. These results can 
be used in conjunction with current and future high-energy observations of bursting neutron stars 
to test general relativity in the strong-field regime. 



I. INTRODUCTION 

Black holes and neutron stars probe the strongest grav- 
itational fields found in the present-day universe. Astro- 
physical systems with properties that are dominated by 
such compact objects are routinely observed throughout 
the electromagnetic spectrum and will soon be detected 
by gravitational wave observatories. It is our hope and 
expectation that such observations will allow for clean 
tests of strong-field general relativity and will demon- 
strate beyond doubt that black-hole candidates are sur- 
rounded by event horizons. 

Testing strong-field gravity and searching for evidence 
of event horizons with current and future observations 
require a theoretical framework within which such ob- 
servations can be understood and possible violations of 
the general relativistic predictions can be quantified. Be- 
cause of the strong gravitational field and the relativis- 
tic velocities found in the vicinities of compact objects, 
the Parametric Post-Newtonian framework [l| , which has 
been very successful in performing weak-field tests, can- 
not be used. 

Current studies are either based on phenomenological 
spacetimes within the context of general relativity Q, 
or employ theories that are derived from parametric ex- 
tensions of the Einstein-Hilbert action 3]. Albeit useful 
in performing null hypothesis experiments and strong- 
field tests, these approaches do not necessarily explore 
the whole range of possible violations of the general rel- 
ativistic predictions. 

In this paper, I show that several observable quanti- 
ties of bursting neutron stars can be calculated without 
requiring the validity of Einstein's field equations. In- 
stead, they can be derived from conservation laws and 
Killing symmetries, assuming only the validity of the Ein- 
stein equivalence principle. This alleviates the need for a 
general parametric theory with which these strong-field 
tests of gravity can be performed. It also allows for a 
direct translation of observable quantities such as red- 
shifts, luminosities, temperatures, and apparent surface 



areas, into values of the metric elements on the surfaces 
of neutron stars (see i)III|) . 

In the calculations reported here, I assume only the va- 
lidity of the Einstein equivalence principle. This is one of 
the fundamental building blocks of gravity theories and 
encompasses the concepts of the equivalence between in- 
ertial and gravitational mass, of local Lorentz invariance, 
and of local position invariance [jj. A violation of the 
equivalence principle has severe implications for many 
non-gravitational experiments. This has allowed solar- 
system and laboratory experiments to place constraints 
on the amplitudes of such violations that reach as low 
as one part in 10 15 [l|, justifying largely the assumption 
that the Einstein equivalence principle is valid even in the 
strong gravitational fields found in the vicinities of com- 
pact objects. Note here that the validity of the Einstein 
equivalence principle does not imply the validity also of 
the strong equivalence principle and hence, in a general 
metric theory, a self gravitating body would not be fol- 
lowing the same geodesies as the test particles considered 
in this study. 

Finally, I assume for simplicity that the neutron stars 
are not rapidly spinning. This allows for the result to be 
expressed in terms of only three parameters of the metric 
of each compact object. The calculation can be extended 
in principle to incorporate effects of rapid rotation. In 
that case, however, the results will depend on at least 
two additional metric elements. 

In fJIIl I describe the formalism used in the calculation 
and justify the relevant assumptions. In mill I calcu- 
late three observables related to the spectra of bursting 
neutron stars that have been used in the past both in 
measuring their masses and radii Q and in constraining 
scalar-tensor gravity theories They are: the gravi- 
tational redshift of atomic lines from their surfaces, the 
touchdown luminosities of radius-expansion bursts, and 
their apparent surface areas during the cooling of the 
thermonuclear bursts. I conclude in QYV\ with a discus- 
sion of the prospect of using the calculations reported 
here in performing strong-field tests of general relativity 
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with X-ray observations of neutron stars. is given by 

1 



' ' j - , —,0,0,0 . (4) 
II. ASSUMPTIONS AND DEFINITIONS obs \TZ ' 



In the following sections, I will calculate a number of 
observable properties of a bursting neutron star, assum- 
ing only the validity of the Einstein equivalence princi- 
ple [lj, without relying on any field equations for the 
underlying gravity theory. I will assume for simplicity 
that the compact object is slowly spinning so that its 
spacetime can be cast in the general form 



ds 2 



-Vfdt 2 + V 2 dr 2 + r 2 dW 



(1) 



I will also assume that the redshift factor, 1Z, and the vol- 
ume correction factor, V, are only functions of the coor- 
dinate radius r with the proper Minkowskian asymptotic 
limits. In this spacetime, I will denote the coordinate 
radius of the stellar surface by r s and the corresponding 
values of the metric elements there by 1Z B and V s . In 
general relativity, as well as in a number of alternative 
gravity theories, 1ZV = 1 in the spacetime external to the 
star and 1Z S V S — 1 in its low-density surface layers. 

Because of the assumption of the validity of the equiv- 
alence principle, matter and photons follow geodesies of 
the spacetime described by the metric (QJ. I will de- 
scribe the motion of a massive particle in terms of its 4- 
velocity = (it* ,u r ,u e ,u^). I will calculate the u* and 
u r components of the 4-velocity using the conservation 
laws that arise from the two Killing vectors with compo- 
nents ^ = (1,0,0,0) and rf = (0,0,0,1) of the space- 
time, i.e., the conservation of energy E = —^u^ = 7?. 2 it* 
and of angular momentum / = r/^u^ = r 2 u^. Without 
loss of generality, I will rotate the coordinate system so 
that the orbit of the particle is at the equatorial plane, 
i.e., I will set sin# = 1 and u e = 0. Finally, I will calcu- 
late the last component of the 4-velocity of the particle 
using the requirement u^u^ = — 1, which leads to 



E 2 



1 / I 2 



1/2 



(2) 



I will describe the motion of a photon in terms of its 
4-momentum — (fc*, k r , k , k^). Because of the same 
Killing vectors and symmetries as in the case of massive 
particles, I can again write fc* = E/1Z 2 , k e = 0, and k^ = 
l/r 2 . Finally, I will calculate the radial component of the 
photon 4- momentum from the requirement k^k^ = 0, 
which gives 



k r = 



( E 2 



\TZ 2 V 2 



V 2 r 2 



1/2 



(3) 



Because the spacetime is asymptotically Minkowskian, 
by assumption, the ratio b = l/E is the asymptotic im- 
pact parameter of the photon at radial infinity. 

In the following, I will also make frequent use of the 4- 
velocity of a static observer at coordinate radius r, which 



III. OBSERVED PROPERTIES OF NEUTRON 
STARS 



Thermally emitting neutron stars are prime candidates 
both for constraining the equation of state of neutron- 
star matter Q and for testing general relativity in the 
strong- field regime [3j. This is especially true for neu- 
tron stars that experience strong thermonuclear flashes 
on their surface layers, which are observed as type I X-ray 
bursts [5|. 

The three observable properties of bursting neutron 
stars that can be used in testing strong-field general rel- 
ativity are the gravitational redshift of a surface atomic 
line, the touchdown luminosity of a radius-expansion 
burst, and the apparent surface area during the cooling 
phases of the bursts. I will now calculate them in detail. 



A. Gravitationally Redshifted Lines 

If the spectrum of a thermally emitting neutron star 
has absorption (or emission) features characteristic of 
atomic transitions, these features will be detected by an 
observer at infinity with a gravitational redshift equal to 



A 



= KI 1 - 1 



(5) 



where Ao is the rest-frame wavelength of the atomic tran- 
sition. For a slowly-spinning neutron star, the emission 
line will be rotationally broadened to a width AA of 



AA ~ 2- 



where Q is the spin frequency of the star. 



(6) 



B. Touchdown Luminosity of a Radius- Expansion 
Burst 



The brightest among the type I X-ray bursts from 
an accreting neutron star show strong spectroscopic evi- 
dence for rapid expansion of the radius of the X-ray pho- 
tosphere [7|. It is widely believed that the luminosities 
of these bursts reach the Eddington critical luminosity at 
which the outward radiation force balances gravity, caus- 
ing the expansion of the surface layers of the neutron star. 
The touchdown luminosities of radius-expansion bursts 
from a given source remain constant between bursts to 
within a few percent, giving empirical verification to the 
theoretical expectation that the emerging luminosity is 
approximately equal to the Eddington critical luminos- 
ity 0. 
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In general relativity, the Eddington critical luminos- 
ity at the surface of bursting neutron star depends on 
its mass and radius. I calculate here this critical lumi- 
nosity for a compact object with an external spacetime 
described by the metric ([J), following refs. 

I define the Eddington limit as the critical flux at which 
the outward radiation force keeps the radial velocity of a 
particle constant, i.e., the one for which 



d 2 ^ 

dr 2 



0, 



(7) 



where r is the proper time for the particle. If / M is the 4- 
force that correspond to the Eddington limit, then from 
the geodesic equation in the radial direction 



f d 2 r 



(8) 



where m is the rest mass of the particle and u M its 4- 
velocity. The only non-zero component of the connection 
that enters the geodesic equation (J8]) is 



K dK 



(9) 



For a particle at rest with a 4- velocity u' 1 = (K~ x , 0, 0, 0), 
the radiation force that corresponds to the Eddington 
limit becomes 



I 



-r = 



1 dK 

KV 2 dr 



(10) 



The outward radiation force is related to the radiation 
flux in the particle's rest frame F r via 



f = oF r 



(11) 



where a is the cross section for interaction. For a particle 
at rest, the radiation flux is given by 



F r = -T {r){t) u t = KT {r){t) , 



(12) 



where T^^' is one element of the energy tensor of the 
radiation field and the parentheses in the indices empha- 
size the fact that no summation over the repeated index 
is implied. The ri-element of the energy tensor of the ra- 
diation field at any coordinate radius r can be calculated 
from (ref. [| eq. (3.32)) 



T®W = 7r/(r) sin 2 a , 



(13) 



where the hats on the indices indicate the fact that these 
are the tetrad components of the tensor. Using the no- 
tation of Ref. [8j] , 



T (t)(r) = J_ T (t)(f) 

KV 



nl(r) sin a 
KV 



(14) 



Here a is the opening angle of the star as viewed by a 
static observer at coordinate radius r, and I(r s ) is the 
frequency-integrated specific intensity of the radiation, 



r 
B 




FIG. 1: The geometry used in deriving the Eddington critical 
luminosity for a slowly spinning neutron star. 



which I consider here to be isotropic within the opening 
angle a. 

In order to calculate the opening angle a, I will follow 
the procedure of ref. as illustrated in Figure 1. Let 
point B be the position of a static observer in the space- 
time with a 4- velocity u^ 1 (eq. Q), and also of two pho- 
tons with 4-momenta u> M and The angle £ between 
the two photons as measured by the static observer is 



cos£ = 1 



(15) 



I will consider one photon to be traveling radially out- 
wards, i.e., with a 4- momentum vector that corresponds 
toZ = (eq-EJ) 



'"-'£,1,0,0 



(16) 



and the other to be is traveling on a non-radial null 
geodesic 



E 



K 2 ' 



1 



K 2 V 2 V 2 r 2 



1/2 



(17) 



where b = l/E. The angle between the two photons as 
measured by the static observer is 



cos 4 = 1 T b 



1/2 



(18) 
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from which it follows that 



leads to 



71 



sin£ 



(19) 



is a constant of motion for each photon. This allows me 
to calculate the opening angle of the star as viewed by 
an observer at coordinate radius r. 

A photon that leaves the stellar surface at point A at 
an angle S with respect to the radial direction is detected 
by the observer at radius r traveling at an angle a with 
respect to the local radial direction. In its motion, the 
photon follows a null geodesic characterized by a single 
impact parameter b and hence 



u rs ■ x r ■ 
o = —— sin o = — sin a 

TZ S TZ 



(20) 



If the stellar surface is at a coordinate radius larger than 
the one of the photon orbit, then the photon with the 
largest angle a detected by the observer is the one that 
emerged from the stellar surface at S = ir/2. If the stellar 
surface is inside the photon orbit then the opening angle 
corresponds to S < ir/2. As a result, the opening angle 
of the source as seen by the static observer is 



sin a < 



2Ln 

n R r 



(21) 



with the equal sign corresponding to stars larger than the 
photon orbit radius. Because I will be looking for small 
deviations from general relativistic neutron stars, all of 
which are larger than the radius of the photon orbit, I 
will consider here only this case. It is interesting that, 
although the details of self-lensing depend on the proper 
volume factor V, the opening angle of the source as seen 
by an observer at infinity depends only on the rcdshift 
factor TZ, as in general relativity. 

The radiation force on a particle becomes 



r = 



aH-Kl{r) ( TZ 



TZV 



(22) 



In order for the particle to remain at a constant radius 
r, this force has to be equal to (see eq. ([TO])) 



1 a 
TZVrn 



7ZirI(r) 



tzZ 



1 dTZ 
TZV 2 dr 



(23) 



Evaluating this expression at the stellar surface, I obtain 



-n 2 7rl(r s 



71 dTZ 
V ~dr 



(24) 



I now connect the specific intensity of the radiation 
field at the stellar radius with that at infinity. Conserva- 
tion of the photon occupation number 



Hrs) 



(25) 



Z(r -> oo) = I(r s )K 4 s , (26) 

from which I can also infer for the luminosity L s — 
4irr 2 nl(r s ) that 

L[r -> oo) = L{r s )n 2 s . 
As a result, condition (|23|) becomes 



(27) 



r oo 

-^E — 



Airmr 2 i 


tz s 


an 


\ 47rm r 2 g c g 


a \ 




dr 


J ' ° (l + *s) 2 



(28) 



where I have defined as the Eddington critical lu- 
minosity as measured by an observer at infinity and the 
effective acceleration at the stellar surface by 



9cS 



1 dTZ 



TZ S V S dr 



(29) 



C. Apparent Surface Area During Burst Cooling 

A typical type I X-ray burst on a neutron star is char- 
acterized by a rapid (~ 1 s) rise and a slower (tens of 
seconds) decay, during which the thermonuclear flash has 
spread throughout the entire neutron star. Observations 
of the cooling tails of multiple type I bursts from a single 
source have shown that the apparent surface area of the 
emitting region, defined as 



S°° = 



4ttD 2 F c , 



CSB 



(30) 



remains approximately constant during each burst and 
between bursts from the same source. In definition (|30p . 
Fc,oo is the measured flux of the source during the cooling 
tail of the burst, T c °° is the measured color temperature 
of the burst spectrum, D is the distance to the source and 
ctsb is the Stefan-Boltzmann constant. In general relativ- 
ity, the apparent surface area during the cooling tails of 
X-ray bursts from a neutron star depends on the stellar 
mass and radius. I calculate here this same quantity for 
an object described by a general spherically symmetric 
metric. 

The color temperature on the surface of the neutron 
star is related to the color temperature measured at in- 
finity by a simple redshift, i.e., 



71* 



(31) 



I can then use equations (fl"2")) , (| 13[) , and (|2"Tj) to calculate 
the radiation flux measured at infinity as 



F Ci00 - TrI(r s )TZ 2 s 



(32) 



obs 



If the spectrum of a neutron star during the cooling 
tail of a burst were a pure blackbody, then nl(r s ) = 



5 



&S~b[T (r s )] 4 . However, the interaction of photons with 
the atmospheric electrons and atoms distort the shape 
of the emerging spectrum. The distortion is typically 
parameterized in terms of a color correction factor 



from its surface, 



/c 



T 

± c 



(33) 



where T c and T c g are the color temperature and the effec- 
tive temperatures of the atmosphere, respectively. The 
value of the color correction factor can be calculated a 
priori and depends on the temperature and the compo- 
sition of the atmosphere as well as on the local accelera- 
tion With this definition, 



7r/(r s ) 



4 

c.s 



./; J 



and the observed ratio -Fc.oo/csb^c 1 ^ becomes 



1 



it \d) 



(Zs + l) 2 



(34) 



(35) 



It is worth emphasizing that the dependence of the ob- 
servable on the coordinate radius of the stellar surface 
and on the redshift is the same as in the case of general 
relativity. 



IV. TESTING GENERAL RELATIVITY WITH 
BURSTING NEUTRON STARS 

In mill I showed that various observables from a slowly 
spinning neutron star depend on three parameters related 
to the coordinate radius of its surface and the values of 
the metric elements there. In general relativity, however, 
all these observables depend on only two parameters: the 
gravitational mass of the compact object and the coor- 
dinate radius of its stellar surface. The difference in the 
number of free parameters between a general metric the- 
ory and general relativity makes a direct test of the lat- 
ter theory possible. Indeed, the three parameters of a 
compact object in a general metric theory can be mea- 
sured (or at least constrained) using high-energy obser- 
vations. If general relativity describes accurately gravi- 
tational phenomena in the strong-field regime, then the 
three parameters cannot be independent but they must 
satisfy a consistency relation, which I calculate below. 
The degree to which consistency can be demonstrated 
will provide a measure of the degree to which general 
relativistic predictions have been tested. 

In general relativity, the two components of the met- 
ric ([1]) at the stellar surface are related to the mass M 
and the radius r s of the compact object by 



1 



2M 



(36) 



Instead of using the mass of the compact object as the 
second independent variable, I will choose the redshift 



Z S; GR 



f 



1 



(37) 



The effective gravitational acceleration in the surface lay- 
ers of a general relativistic star then becomes 



Scff.GR 



1 

2^7 



z s (z s + 2) 
z s + l 



(38) 



I can now define a parameter r\ to measure deviations 
from general relativity by 



1 

2r s 



z s (z s + 2) 



Z S + 1 



(39) 



This parameter quantifies the degree to which the grav- 
itational acceleration on the surface of a star can be de- 
termined entirely by its mass and radius, as predicted by 
general relativity. 

With the above definitions, I can write the three ob- 
servable quantities for a bursting neutron star of known 
distance, D, as 



SX 



D 



47rm p r s 
(1 + X)a T 

V2 „ 



fi 



z s (z s + 2) 

(i + z s r 

:(z s + l) , 



(40) 
V (41) 



(42) 



where I have introduced the apparent source radius i? a pp- 
I have also assumed that the interaction between the elec- 
trons and the photons is due to coherent Thomson scat- 
tering and denoted the hydrogen mass fraction of the 
accreted material by X. 

It is important to note that, for a general metric theory 
of gravity, the first and the third observables depend on 
the coordinate radius of the neutron-star surface, r B , and 
on its gravitational redshift, z s , in the exact same way as 
in general relativity. On the other hand, the Eddington 
luminosity depends also on the parameter rj, which mea- 
sures possible deviations from general relativity. There- 
fore, it is only through measurements of the touchdown 
luminosity of radius expansion bursts that general rela- 
tivity can be tested with bursting neutron stars. 

In general relativity, t^qr = 1> by definition. On the 
other hand, in other gravity theories, this parameter de- 
pends also on the additional gravitational degrees of free- 
dom. For example, in a scalar-tensor gravity as the one 
described in ref. Q), rj depends on the magnitude of the 
scalar field at the surface of the neutron star, which in 
turn depends on the mass distribution inside the star and 
its coupling to the scalar field. 

Not all combinations of values of the three observables 
are possible in general relativity. Buchdahl's theorem 
states that no spherical, general relativistic star can ex- 
ist with a mass-to-radius ratio larger than M/R > 4/9. 
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FIG. 2: Solid lines: Contours of constant Eddington luminosity 
at infinity, in units of 10 38 erg s — 1 , as a function of the coordinate 
radius and the redshift of the surface of the bursting neutron star 
Dashed lines: Contours of constant apparent radius of the neutron 
star, in units of 1 km, as inferred from the cooling tail of a burst, 
on the same parameter space. In this calculation the atmosphere 
of the neutron star was assumed to consist of pure hydrogen and 
the color correction factor was set to 1.34. 



This corresponds to a maximum possible value of the 
redshift from the surface of a general relativistic star of 
z s < 2. The apparent radius during the cooling tails 
of bursts depends monotonically on the redshift and, 
therefore, reaches the highest possible value in general 
relativity at z s = 2. Finally, the Eddington luminos- 
ity is a non-monotonic function of redshift and peaks at 
z s = V3 - 1 ~ 0.73. The maximum possible values of 
the three observables for a general relativistic star are, 
therefore, 



SX 



< 2 



GR 



T OO 

L E,GR 



< (V3-1) 



47rmr s 
(l + X)<r T 



R 



app,GR ^ n 2 
7c 



(43) 
(44) 
(45) 



Violation of any of the above inequalities will also sig- 
nify that general relativity does not describe accurately 
gravitational phenomena in the strong-field regime. 

A second test of general relativity can be obtained us- 
ing a slowly-spinning, bursting neutron star at a known 
distance that accretes matter of known composition and 
shows radius-expansion type I X-ray bursts. Assuming 
the validity of general relativity, equations (|4Tj) and (142)) 
can be inverted to yield the coordinate radius of the 
neutron-star surface and its gravitational redshift (see 
Fig. [2]). However, because the Eddington luminosity does 
not have a monotonic dependence on redshift, this in- 
version procedure is not always possible. For example, 
Figure shows that an apparent radius of 5 km and an 
Eddington luminosity of 2 x 10 38 erg s _1 are impossible 



in general relativity. 

In order to calculate the range of simultaneously de- 
termined values of the Eddington luminosity and of the 
apparent radius that are consistent with general relativ- 
ity, I solve equation (|4ip for r s and substitute it into 
equation (|42|1 . After a small rearrangement of terms, the 
result is 



(l + A> T Lg GR _ z s (z s + 2) 
47rm p / c 2 i? app!GR (z s + l) 4 



(46) 



The right-hand-side of this last equation has a maximum 
value of 1 /4 at z s = \/2 — 1 , so that 



L 



E.GR 



< 



{l + X)a T 



R 



app,GR 



which I can rewrite in CGS units as 



L™ GR < 1.9x10 



:i,S 



fc 



1.34 



R 



app,GR 



(47) 



erg s 



1 + X J V 10 6 cm 

(48) 

Violation of inequality (j48|) will also be evidence for new 
gravitational physics, not described by general relativity. 

The optimal test of general relativity will occur from 
the detection of gravitationally rcdshiftcd lines from a 
neutron star that is slowly spinning, lies at a known dis- 
tance, and shows radius-expansion type I X-ray bursts. 
In this case, equations (|40|) - (|42|) can be combined to give 



(1 + AW / F c . 



mfc 



\ u c,c 



-1/2 



(1 + SX/X ) 4 
+ SX/\ )(SX/\ 



-D 



(49) 

All the quantities on the right-hand side of this last equa- 
tion are either observable or can be calculated from mod- 
els of neutron-star atmospheres. The parameter r\ is a 
measure of the degree to which general relativity accu- 
rately describes the strong gravitational fields found in 
the vicinities of neutron stars. 

A successful application of the line of arguments dis- 
cussed here relies on the detection of radius-expansion 
bursts and of gravitationally redshifted lines from sources 
with known distances. Radius expansions bursts have 
been detected to date from at least 31 bursting neutron 
stars [l2j], several of which lie in globular clusters and, 
therefore, have also well constrained distances The 
touchdown luminosities and apparent radii can be in- 
ferred from current data to an accuracy of ~ 5% that 
is limited by systematic uncertainties only Q- The 
distances to the globular clusters that contain known 
bursters are known to an accuracy of ~ 5 — 20%. As- 
suming the latter accuracy to be 10%, and assuming also 
that the composition of the accreting material can be es- 
timated to the same accuracy by long-wavelength obser- 
vations of the accretion flow, the proposed test can lead 
to a measurement of the parameter r\ to within ~ 15%. 

Note, however, that gravitationally rcdshiftcd atomic 
lines have been reported only for one, slowly rotating, 
burster, EXO 0748—676. Future observations of burst- 
ing neutron stars with telescopes that combine a high 
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spectral resolution and a large collecting area, such as tests of general relativity in the strong-field regime with 
Constellation-X and XEUS, will allow for quantitative bursting neutron stars. 
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